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We use a semiclassical three-fluid thermodynamic model to evaluate the phenomena of spatial 
demixing in mesoscopic clouds of fermionic and bosonic atoms at high dilution under harmonic 
confinement, assuming repulsive boson-boson and boson-fermion interactions and including account 
of a bosonic thermal cloud at finite temperature T. The finite system size allows three different 
regimes for the equilibrium density profiles at T — 0: a fully mixed state, a partially mixed state 
in which the overlap between the boson and fermion clouds is decreasing, and a fully demixed state 
where the two clouds have zero overlap. We propose simple analytical rules for the two cross-overs 
between the three regimes as functions of the physical system parameters and support these rules by 
extensive numerical calculations. A universal "phase diagram" expressed in terms of simple scaling 
parameters is shown to be valid for the transition to the regime of full demixing, inside which we 
identify several exotic configurations for the two phase-separated clouds in addition to simple ones 
consisting of a core of bosons enveloped by fermions and vice versa. With increasing temperature 
the main role of the growing thermal cloud of bosons is to transform some exotic configurations 
into more symmetric ones, until demixing is ultimately lost. For very high values of boson-fermion 
repulsive coupling we also report demixing between the fermions and the thermally excited bosons. 
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I. INTRODUCTION 

After the realization of Bose-Einstein condensation (BEC) in atomic gases one of the most challenging endeav- 
ours in experiments on cold atoms is the cooling of fermionic isotopes of alkali atoms down to expected the superfluid 
transition. Theoretical estimates based on the BCS model Q indicate that the temperature for the superfluid transi- 
tion should be much lower than the Fermi temperature Tf, which is of the order of the BEC transition temperature. 
In fact, cooling fermions is harder than cooling bosons: the main difficulty arises from Fermi statistics, as s-wave 
collisions between spin-polarized fermions in a magnetic trap are forbidden by the Pauli principle. A common strategy 
uses sympathetic cooling, which is based on s-wave collisions between the fermions and a second gaseous component 
made either of fermions in a different internal state or of bosons. The latter choice seems to minimize the effects of 
Pauli blocking, which limit the process of cooling two fermionic components ||. Several experiments are currently 
in progress on trapping and cooling various boson-fermion mixtures, i.e. 6 Li- 7 Li @,|], 6 Li- 23 Na § and 40 K- 87 Rb @. 
The lowest temperature attained in these experiments so far is about 0.2 Tf. 

From the theoretical point of view a mixture of condensed bosons and fermions in the normal state is already an 
interesting system to study, because it can show spatial demixing of the two components Q . In the homogeneous gas 
at T = this is an example of a quantum phase transition, that is a phase transition induced by the interactions [0]: on 
increasing the boson-fermion repulsion the system minimizes its total energy by placing the bosons and the fermions 
(or boson-fermion mixtures of different compositions) in different regions of space Jl0| , even though this implies a high 
cost of kinetic energy at the interface. Boson-fermion demixing in a quasi-spherical trap has been studied by Nygaard 
and M0lmer pip, while the conditions for demixing inside a spherical trap have been set out by two of us E^] for 
. £h ! N\, — Nf and by Miyakawa et al. [|l3| for Nf <C iV&, with Nb and Nf the numbers of bosons and fermions in the trap. 
The conditions for observing demixing in the Paris experiment on a 6 Li- 7 Li mixture in a cigar-shaped trap § has 
been analyzed by Akdeniz et al. JTif . There are, of course, important consequences of finite system size on demixing 
in trapped mixtures, relative to the case of a homogeneous one. In particular, the transition is spread out in a finite 
system and the anisotropy of the trap acts differently on the two types of atoms (expressions such as phase transition 
and phase diagram will nevertheless continue to be used for brevity in the following). Locating the onset of partial 
demixing is also relevant to the practicalities of fermion cooling, since at that point the diminishing overlap between 
the two clouds will start reducing the effectiveness of the collisional transfer. 

In this paper we give the general conditions under which phase separation of bosons and fermions occurs inside 
a harmonic trap at varying anisotropy and relative numbers of the two components. Within a mean-field model we 
make an extensive study of the equilibrium density profiles of the bosonic and fermionic component, ranging from the 
fully mixed state at small values of the scattering lengths to the partially demixed state at intermediate values of the 
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coupling constants and finally to the regime of full phase separation. We show that the results can be understood on the 
basis of scaling laws and simply summarized into a phase diagram, which is expressed in terms of scaling parameters 
of the system at zero temperature. In the phase-separated regime we find several configurations having different 
symmetry and topology, some of which are metastable, and investigate the role of the anisotropy of the confinement 
in determining the minimum-energy configurations. We then extend our analysis to finite temperature, finding that 
on increasing temperature some phase-separated configurations turn or decay into others of higher symmetry before 
phase separation disappears. Finally, we show that phase separation between the fermions and a bosonic thermal 
cloud is also possible in principle on further increase of the boson-fermion scattering length. 

The paper is organized as follows. In Sec. || we describe the model that we have used and give the limits of its 
applicability, together with a schematic description of the numerical method employed to find the phase-separated 
configurations. Sec. Ill summarizes the conditions under which phase separation occurs at T — 0, illustrates various 



configurations which can be found in the phase-separated regime and gives the phase diagram for the lowest energy 
configurations. Sec. IV illustrates the effect of temperature on the configurations obtained for both small and large 
boson-fermion scattering length. Finally, Sec. [v| gives a summary and some concluding remarks. 



II. THE METHOD 

We describe the boson-fermion mixture by means of the particle density profiles, which are n c (r) for the condensate, 
n nc (r) for the bosonic thermal cloud and nf(r) for the fermions. The components are subject to axially symmetric 
confining potentials given by 



(1) 



where my are the atomic masses, u>bj the trap frequencies and X b j the trap anisotropies. We evaluate the density 
profiles in a mean-field model using the Thomas-Fermi approximation for the condensate and the Hartree-Fock 
approximation for the other clouds |l5| , |l6|] . 

The Thomas-Fermi approximation assumes that the number of condensed bosons is large enough that the kinetic 
energy term in the Gross-Pitaevskii equation may be neglected |0| . It yields 



,(r) = \l* b -V b e * t (T) 



2gn nc (r) - fn f (r)]/g 



(2) 



for positive values of the function in the square brackets and zero otherwise. Here, [ib is the chemical potential of the 
bosons and the coupling constants are g = 4irh 2 abb/rrib and / = 2ith 2 a,bf /m r , with ay, and a^f the boson-boson and 
boson-fermion scattering lengths and m r = rribmf / (nib + nif) the reduced mass. The Hartree-Fock approximation, 
on the other hand, treats the thermal boson cloud and the fermion cloud as ideal gases subject to effective potentials 
V b e j f , that is 
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and 



V b eff (r) 



V b ext (r) + 2gn c {v) + 2gn nc (r) + fn f (v) 



Vf f {v) = Vf xt {v) + fn c (r) + fn nc (r). 
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The chemical potentials fi^j are determined by requiring that the volume integrals of the densities n c (r) +n nc (r) and 
n/(r) should be equal to the numbers Nb and Nj of particles. 

The model is valid for the bosons when the diluteness condition n c a\ b -C 1 holds and the temperature of the 
mixture is sufficiently below the condensation temperature. The fermionic component has been taken as dilute spin- 
polarized gas, for which the fermion- fermion s-wave scattering processes are inhibited by the Pauli principle and 
p-wave scattering is negligible fl8|| . The condition kfdbf <C 1 with kf the Fermi wave number should hold in the 
mixed regime, but this is not a constraint in the regime of phase separation where the boson-fermion overlap is rapidly 
dropping. 

The density profiles at fixed numbers Nb and Nf are determined numerically by a self-consistent solution of the above 
equations. The profiles are obtained in two self-consistency loops: first, Eqs. (g) and (||) are solved at fixed chemical 
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potentials starting from some initial guesses for the densities; then the chemical potentials are found iteratively by a 
standard algorithm for multi-variable minimization. Whereas in the mixed state the density profiles at convergence do 
not depend on the initial choice, different configurations are found in the phase-separation regime for the same values 
of the parameters by varying the initial conditi ons. We have used this fact to search for several possible metastablc 
configurations, which will be illustrated in Sec. [II B below. 



III. PHASE DIAGRAM AT ZERO TEMPERATURE 



As is evident from the preceding section, a boson-fermion mixture under confinement is characterized by a large 
number of parameters. The aim of this section is to analyze what is the effect of varying them independently and to 
provide a unified understanding of the results using scaling laws. 

The interaction energy grows on increasing the values of the boson-boson and boson-fermion coupling constants. 
There exists a threshold beyond which the total energy is minimized by a configuration in which the two components 
are spatially separated. The transition to the phase-separated regime is smooth owing to the finite size of the confined 
gas, and three different states are indeed recognizable: mixed, partially demixed, and fully demixed. Depending on 
the value of the relative strength of the two coupling constants, it may happen that either the fermions or the bosons 
are pushed away from the center of the trap. 



A. Partial demixing 



The boson-fermion interaction energy is proportional to the overlap between bosonic and fermionic clouds and at 
T = we have 

E mt = f J d 3 rn c (r)n f (r). (6) 

We locate the onset of partial demixing by looking at the point where the boson-fermion overlap starts to decrease 
with increasing abf/abb- 

In the following we focus on the simplest case where only one scattering length is varied, which also seems realistic 
from an exper imental point of view. Then the most efficient way to approach phase separation is by changing a^f 
(see Sec. Ill B below). The behavior of the interaction energy as a function of the ratio abf/abb at fixed abb is shown 



in Fig. 1(a) 



The position of the maximum of the interaction energy at fixed abb can be estimated analytically from the 
condition dEi nt /df — 0, by using the approximate expressions nbj ~ Nbj /(4ttR 3 y/3) and the values for the 

cloud radii Rf and Rb as obtained in the absence of boson-fermion interactions, Rf — (48iV//A/) 1 / 6 a/ and 

Rb = {lhXbNbabb/ 'ab) 1 ^' at>/ 'A^ 3 with abj — (h/mbjLUbj) 1 ^ 2 . The expression for the position of the maximum is 
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We recognise in Eq. (Q) a geometric combination of two scaling parameters: c\Ny 2 / iV^ is dominant in the case 



/ 

/N^ 3 was previously identified |l^] in the regime Nf <C iVj,. The predictions obtained from 



r 2 /5//\ r i/3 



N b < Nf while c 2 N b 

Eq. (0) are indicated in Fig. 1(a) by vertical arrows. There clearly is very good agreement between the analytical 
estimate and the numerical results. 
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B. Full demixing 



The regim e of full demixing is taken to be reached when the boson-fermion overlap becomes negligible. As is shown 



in Fig. 1(b) , within the three-fluid model the transition point does not depend on the number of bosons in the trap. 
This can also be predicted from the results obtained by Viverit et al. [ [To| for the homogeneous mixture, by using the 
valu es of the densities taken at the center of the trap. Within the same approximation for the fermion density as in 



Sec. [II A| we obtain the condition for phase separation at T = as 

akfa bb > I — ) , (10) 

\ a bfJ 

where k f = (48A/A 'f) 1 / 6 / 'a f and a = [3 1 / 3 /(2ir) 2 / 3 ]m b m f /(4m?,). Equation © is valid for different numbers of 
bosons and fermions as well as for different atomic masses and trap frequencies, and agrees very well with the phase- 
separation criterion given earlier by two of us [l2| for the case where mj = nib and u>f = ui b - It confirms that the 
two relevant parameters for describing the transition to the fully demixed regime are abf/abb and akfabb- We also 

1/12 

immediately find that the value of abf at the transition point scales as N* at fixed abb- this result is illustrated 



in the inset in Fig. 1(a) 



Equation ([L0|) has been verified by the numerical solution of the three-fluid model for a variety of different sets 
of parameters. First of all we have studied the location of phase separation in parameter space by varying the 
scattering lengths over much wider ranges than allowed by the usual off- resonant values (e.g. in the 6 Li- 7 Li mixture), 
as might be attained experimentally by exploiting optically or magnetically induced Feschbach resonances p0[ . For 
these calculations we have chosen the geometrical parameters of the traps and the numbers of atoms as in the Paris 
experiment on 6 Li- 7 Li g: namely N b ~ N f = 10 4 , w fc /27r = 4000 Hz, uj f /2tt = 3520 Hz and A 6 ~ A/ ~ 1/60. The 
results are summarized in a quantum phase diagram at T = in the plane {abf/abb, kfabb}, which is reported in 
Fig. |^. This contains all the lowest-energy configurations that we have found, with the inset showing an enlargement 
of the region near the origin of the above-mentioned plane. The curve drawn inside the inset in Fig. |^ is obtained 
from Eq. (|lo| ) and is in full agreement with the numerical results in separating the region of full demixing from the 
mixing or partial-demixing regions. 

We turn to a detailed account of the configurations corresponding to the various symbols and letters in the phase 
diagram in Fig. ^. Below the phase separation line in the inset two different types of configurations are found, one 
with a core of bosons in partial overlap with an envelope of mainly fermions (triangles down) and a specular one with 
a core of fermions enveloped by mainly bosons (triangle up). An interesting configuration that we have observed in 
the latter case shows coexistence of a mixed state with a purely bosonic phase (see Fig. ||): this combination is not 
stable in the homogeneous mixture |l0| . 

In the regime of phase separation a variety of different configurations are observed, including several that require 
a break of the symmetry imposed by the confining harmonic potentials. The letters a-d in the phase diagram in 
Fig. H give the locations of four such energetically stable structures, which are shown in Fig. [| together with some 
other metastable structures. The symmetric configuration with boson inside and fermion outside, indicated by (a), 
is obtained in Fig. |^ with the choice abb = 600 Bohr radii and abf — 10 abb, while the complementary configuration 
(b) formed from fermions inside and bosons outside is obtained with a b b = 64600 Bohr radii and abf — abb/2. 
Other configurations obtained with the same scattering lengths as the symmetric structure of type (a) are shown in 
Fig. ^: these are (c) a boson torus inside a fermion cloud, (d) a threefold symmetric structure formed from fermions 
surrounded by a shell of bosons inside a fermion envelope, (e) an asymmetric structure in which the bosons are shifted 
away from the center of the trap along the z axis, and (f) a "sandwich" formed by bosons inside a fermion cloud. 
Finally the structure indicated by (g) in Fig. which consists of a torus of fermions around an elongated boson core, 
is found with the choice a bb = 6250 Bohr radii and a&/ = 2 abb- 

With the exception of (a) and (b) the configurations shown in Fig. [|are very different from the density distributions 
that are met in the mixed phase, and thus could be identified in an experiment by looking at the column-density 
pictures of the atomic clouds. As an illustration we give in Fig. the column-density images corresponding to the 
case (c) of a bosonic torus, and in Fig. [j] those for the case (f) of a boson sandwich. 

Finally, we have examined the role of the anisotropy of the confinement on the configurations of lowest energy in 
the phase-separated regime. To this end we have evaluated the stability of the lowest-energy configuration by varying 
the anisotropy from A/ = l/60toA/ = l while keeping the ratio Xb/Xf — 1, with the choice abb — 600 Bohr radii and 
abf — lOa^f, for the scattering lengths. The most stable configuration remains the symmetric structure (a) except 
that for Xf ~ 1 the structure (d) has the lowest energy. For this particular choice of parameters the relative difference 
in energy between the metastable configurations does not vary much with A / . 
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IV. THE ROLE OF TEMPERATURE 



At finite temperature an increasing fraction of the bosons populates the thermal cloud, thus depleting the conden- 
sate. In harmonic confinement the thermal cloud is wider and more dilute than the condensate, and thus is much 
less interacting with the fermions. As a result the thermal cloud is not yet in phase separation when the condensate 
already is. 

At low temperatures the non-condensate fraction is small and its presence does not affect the description of the 
phase separation between the condensate and the fermions that we have given in the preceding section. In particular 
we have verified that at T of the order of 0.2 Tf, as are reached in current experiments, all the exotic configurations 
(c) to (g) remain possible. They are not found at higher temperatures: we show as an example in Fig. ^ the behavior 
of the boson torus configuration (c) as a function of temperature. For 0.2 Tf < T < 0.3 Tf a part of the condensed 
bosons moves from the torus to the center of the trap and for T > 0.3T/ the phase separation is lost. 

At still higher temperatures the thermal cloud starts to be consistently populated and we have found that phase 
separation between the fermions and the entire bosonic cloud becomes possible for very large values of dbf/o-bb [ p"4[ . 
In this case, as is illustrated in Fig. g, the fermion cloud exhibits a central hole having the size of the bosonic thermal 
cloud. 



V. SUMMARY AND CONCLUDING REMARKS 

We have performed an extensive analysis of the phenomenon of phase separation for a boson-fermion mixture under 
confinement on the basis of the atomic density profiles of the two components. First of all we have pointed out that 
three regimes can be identified owing to the finite size of the system: a mixed phase where the bosonic and fermionic 
clouds overlap, a partially demixed state where the overlap is decreasing, and a phase-separated regime where there is 
no overlap between the two components at zero temperature. We have given analytical expressions for the positions 
of the two cross-overs between these regimes in terms of the physical parameters of the system, and verified these 
conditions by an extensive numerical study based on varying the coupling constants, the numbers of atoms, and 
the geometry of the trap. General expressions for the transition lines are particularly useful in view of the several 
experiments in progress [§-0, which are run with different experimental parameters such as the relative numbers of 
bosons and fermions, the atomic masses, and the trap anisotropies. 

The transition to the partially phase-separated regime, that we have located where the boson-fermion interaction 
energy starts to decrease on increasing the boson-fermion scattering, is important for the experiments on fermion 
cooling: the boson-fermion collision rate is decreased by a geometrical factor which is proportional to the overlap 
between the two clouds. At fixed boson-boson scattering length we have given a simple model for the transition in 
terms of a combination of two scaling parameters, which characterize the cases Nf, <C Nf and Nf <C A^. 

The transition to full demixing has been illustrated by a universal phase diagram expressed in terms of scaling 
parameters of the system. It is remarkable that in our model the location of this transition does not depend on 
the number of bosons. In the region of phase separation we have found several configurations with various topology 
and identified the minimum-energy configurations. We have also investigated the effect of increasing temperature on 
phase separation. Some exotic configurations that arc found in the phase diagram at T = turn into others of higher 
symmetry before the condensate-fermion demixing is lost. 

The main limitations of the present model come from the two basic approximations that we have employed: the 
Thomas- Fermi approximation for the condensate and the mean-field Hartree-Fock approximation for the other gaseous 
clouds. The Thomas- Fermi approximation breaks down for a very small condensate or for negative values of the boson- 
boson scattering length, and one would need to solve the full Gross-Pitaevskii equation. A description beyond mean 
field is required when the dilution parameters increase and for a trapped gas this is still an open problem. The 
analogue of the Bcliacv expansion at T — has recently been derived for the homogeneous boson-fermion mixture 
@. 

Of course, the evolution of the mixture towards phase separation modifies not only its equilibrium properties, but 
also its dynamical properties and collective excitation spectrum. The approach to the transition will be signalled 
by a softening of the frequencies of surface or bulk modes having the appropriate symmetry for a phase-separated 
configuration to be attained. Experiments measuring dynamical properties could therefore provide an alternative 
method for revealing the onset of demixing. Calculations on these spectra are in progress. 
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(a) (b) 

FIG. 1. Boson-fermion interaction energy at T — (in units of NhuJb where N = Nf + Nb) as a function of the adimensional 
ratio abf/dbb for a^b = 5.1 ao with ao being the Bohr radius. The other parameters are uob/2ii = 4000 Hz, Uf /2tt — 3520 Hz, 
At ~ Xf ~ 1/60, and atomic masses corresponding to the 6 Li- 7 Li mixture, (a) For N = 2.1 • 10 5 , with iV/,1 = 10 4 (continuous 
line), = 1-05 ■ 10 J (dotted line) and = 2 • 10 5 (dotted-dashed line). The arrows indicate the positions of the maxima 
according to Eq. (Q). In the inset are shown the same curves in the same vertical units, after rescaling the second and third 
curve along the horizontal axis by factors (N fi2 /N fA ) 1/12 and (N f/A /N fA ) 1/12 . (b) For fixed number of fermions N f = 10 4 and 
different numbers of bosons Nb — 2 • 10 5 (continuous line) and iVj, = 10 4 (dashed line). 
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FIG. 2. Phase diagram at T — 0, in the plane defined by the adimensional parameters a^fc/ and abf /abb- The bottom left 
corner of the figure is shown enlarged in the inset. The continuous line in the inset corresponds to the condition of phase 
separation given by Eq. (|lo|). In the mixed-phase regime we have the following configurations: A fermions at the center: V 
bosons at the center; T bosons an the center in an almost demixed configuration. For the phase-separated regime we have 
used the following notations (see Fig. ^): (a) symmetric configuration with fermions outside; (b) symmetric configuration with 
bosons outside; (c) boson torus inside a fermion cloud; (d) three-component symmetric configuration with fermions outside and 
in the central core. All the configurations in this phase diagram have been evaluated at zero temperature, using Nb = Nf = 10 4 , 
Ubjl-K = 4000 Hz, uJf/2-K — 3520 Hz, A;, ~ Xf ~ 1/60, and atomic masses corresponding to the 6 Li- 7 Li mixture. 
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FIG. 3. Section of the density profiles (in units of a// ) as a function of z/ab for fermions (solid line) and bosons (dashed line) 
in a configuration where a mixed phase and a purely bosonic one coexist (A in Fig. ^|). This configuration has been obtained 
at T = with abb = 323000 ao and atf = 32300 ao. The numbers of atoms and all other parameters are as in Fig. |^. 
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FIG. 4. Configurations in the phase-separated regime for Nf = Nb = 10 4 particles at T = 0. (a) symmetric configuration 
with fermions outside; (b) symmetric configuration with bosons outside; (c) boson torus inside fermion cloud; (d) threefold 
symmetric configuration; (e) asymmetric configuration; (f) boson sandwich inside a fermion cloud; and (g) torus of fermions 
around an elongated core of bosons. The structures (a) and (c-f) are for the same values of abb ~ 600 ao and abj — 6000 ao, 
and their energies per particle are 4.24, 4.71, 5.38, 5.50, and 5.85 tkiJb- The configuration (b) was obtained for abb = 64600 ao 
and abf = 32300 ao and the value of its energy per particle is 13.25 hujb- The configuration (g) was obtained for a^ = 6250 ao 
and abj = 12500 ao and the value of its energy per particle is 9.27 hub- All other parameters are as in Fig. ^. 
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FIG. 5. Column densities (in units of a^ 2 ) in the radial plane (top) and in the {£, z| plane (bottom) for the fermionic 
cloud (left) and the bosonic cloud (right) of the bosonic torus (configuration (c) in Fig. H, obtained with abb = 600 ao and 
dbf = 6000 ao). The size of the figures in the {x,y} plane and in the {x,z} plane is llfim x 11/im and llfim x AOOfim, 
respectively. The other parameters are as in Fig. El 




FIG. 6. Column densities in the {x, z} plane (in units of a^ 2 ) for the fermions (left) and for the condensate (right) in the 
sandwich configuration (f) in Fig. ^, obtained with abb ~ 600 ao and abf = 6000 ao. The size of the figures is llfim x 400/im. 
The other parameters are as in Fig. bl 
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FIG. 7. Sections of the density profiles at different temperatures (in units of a^~ 3 ) for the fermions (solid line) and for the 
condensate (dashed line) as functions of the radial coordinate r/ab for the configuration (c) in Fig. ^j, obtained with abb = 600 ao 
and dbf = 6000 ao. The values of the other parameters are as in Fig. 0. 
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FIG. 8. Column densities in the {x,y} plane (in units of a^ 2 ) for fermions (left) and bosons (right) in the fully demixed 
regime at T = 0.6 Tf, for abb ~ 5.1 ao and a,bf = 2- 10 s ao. The values of the other parameters are as in Fig. |^. In the right panel 
the thermal bosons are not visible because of the high density of the central condensate. The size of the figures is 11/im x 11/im. 
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